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We analyze the physical mechanisms leading either to syn- 
chronization or to the formation of spatio-temporal patterns 
in a lattice model of pulse-coupled oscillators. In order to 
make the system tractable from a mathematical point of view 
we study a one-dimensional ring with unidirectional coupling. 
In such a situation, exact results concerning the stability of 
the fixed of the dynamic evolution of the lattice can be ob- 
tained. Furthermore, we show that this stability is the re- 
sponsible for the different behaviors. 
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I. INTRODUCTION 

Among the collective phenomena that are currently at- 
tracting the interest of the scientific community one of 
the most relevant concerns the synchronization of the 
temporal activity of populations of interacting nonlin- 
ear oscillators, due to its ubiquity in many different 
fields of science. Experimental evidences of this phe- 
nomenon have been reported for centuries [EJ but in the 
last decades the advance in the comprehension of its na- 
ture has allowed to develop a theoretical description. 
In this context, several successful ideas has been sug- 
gested. An interesting approach proposed in Q-Q] has 
been shown to be useful to describe the dynamic evo- 
lution of the population. The idea consist in modeling 
the system as an assembly of phase oscillators interact- 
ing through continuous-time couplings. For sufficiently 
large coupling strength the system may undergo a phase 
transition from incoherence to spontaneous mutual syn- 
chronization. More challenging from a theoretical and 
realistic point of view is to consider networks of pulse- 
coupled oscillators which may account for the behavior 
of heart pacemaker cells, integrate and fire neurons, and 
other systems made of excitable units. The intrinsic non- 
linearities associated to these models make their dynam- 
ical evolution more difficult to describe and only in the 
last years real advances have occurred |^-^| . 

Up to now, almost all the theoretical approaches have 
been centered around mean-field models or populations 
of just a few oscillators. From these studies it is possible 
to investigate the mechanisms relevant for the formation 
of assemblies of synchronized elements as well as other 



spatio-temporal structures. However, these mean-field 
descriptions are, in many cases, far from reality and other 
methods where the specific topology or geometry of the 
system, as well as the precise connectivity between units, 
must be considered because their effects may be crucial. 
In a such new world many points remain open. In par- 
ticular the majority of works rely on simulations showing 
the outstanding richness that a low dimensional system 
of pulse-coupled oscillators may display. Some examples 
are self-organized criticality, chaos, quasiperiodicity, etc. 

Unfortunately, a rigorous mathematical description of 
these systems is still missing. Some of the theoretical 
papers appeared in the scientific literature proof the sta- 
bility of some behaviors O,0 but they do not explore 
the mechanisms leading to them. The goal of this paper 
goes in both directions: the analysis of the mechanisms 
which are responsible for synchronization and formation 
of spatio-temporal structures, and, as a complement, to 
proof under which conditions they are stable solutions of 
the dynamical equations. Since our motivation is to ana- 
lyze the essence of the problem we have considered a ID 
model which will allow us to illustrate the ideas in a very 
clear way. In spite of this apparent simplicity, this sys- 
tem displays a rich set of behaviors, that depend on the 
specific values of the parameters of the model, which has 
been observed in lattices with higher coordination num- 
bers. 1 12 1^] Notice that populations of ID pulse coupled 
oscillators are currently of great interest in some areas of 
science. As an example let us mention that for a certain 
type of cardiac arrhythmia there is an abnormally rapid 
heartbeat whose period is set by the time that an exci- 
tation takes to travel a circuit. This observation can be 
explained by modeling appropriately the circulation of a 
wave of excitation on a one-dimensional ring fl4|| . In a 
different context, synchronization and periodic states of 
ID populations of phase-locked loops have been recently 
investigated, (jll 15 

The structure of this paper is as follows. In Sec. II we 
describe the system as well as the notation used through- 
out the paper. Sections III and IV are devoted to analyze 
the simplest cases of three and four oscillators, respec- 
tively. In Sec. V we study the general case, whereas in 
the last section we present our conclusions. 
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II. THE MODEL 

Let us consider a system formed by a population of 
(N + 1) oscillators distributed on a ring. The state of 
each oscillator is described by its phase which increases 
linearly in time, until one of them reaches a threshold 
value that, without loss of generality, we have considered 
equal to <fith = 1- When this happens the oscillator fires 
and changes the state of its rightmost neighbor according 
to 



> 1 => 



^ ^ 

*>i+i — > 4>i+i + £<A;+i = M^i+i 



Vi = 0,...,JV (1) 

subjected to periodic boundary conditions, i.e. N+l = 0, 
and where e denotes the strength of the coupling. From 
an effective point of view, the pulse-interaction between 
oscillators, as well as the state of each unit of the sys- 
tem, can be described in terms of changes in the phase, 
or in other words, in terms of the so called phase re- 
sponse curve (PRC), fi<p in our case. Behind this fact 
one assumes that the phase shift elicited by an impulse 
affects the period of a given unit in the current time in- 
terval but not in future intervals. In this paper we have 
also considered a linear PRC |16|. In practise, however, 



this condition can be relaxed since a nonlinear PRC does 
not change the qualitative behavior of the model pro- 
vided the number of fixed points of the dynamics is not 
altered. Moreover, a linear PRC has the advantage of 
making the system tractable from an analytical point of 
view. 

Let us describe the notation used in the paper. The 
population is ordered according to the following criterion: 
The oscillator which fires will be always labeled as unit 
and the rest of the population will be ordered from 
this unit clockwise. After the firing, the system is driven 
until another oscillator reaches the threshold. Then, we 
relabel the units such that the oscillator at <f> = 1 is now 
unit number 0, and so on. The whole process can be 
described through a suitable transformation. This fact 
will enable us to study the origin of different structures 
in a very simple way. Our strategy has been to trace 
the phases of the oscillators after each firing and then 
to construct return maps either of a complete cycle, in 
which all the oscillators fire exactly once, or after a single 
firing + driving process (FD). Let us clarify this point 
mathematically. The first step is to construct the matrix 
of the transformation for a FD. To illustrate the situation 
let us consider the general transformation for a 'jump' 
n between two successive firings, distinguishing between 
n = 1 





H<t>i 



1 - nfc = 

1-^01+02 =(j) 



4>N 

and n > 1 



h = 1 
h 



4>N 



o 



i - fi<f>i + . 



1 - fi<j>i +4>N = 



fn+1 



fn+1 



1 - 

1 - 

1 - 

1 - 



4>2 



N-l 



J N-n+l 
V-ri+2 
V-ri+3 



J N-n+i+l 



fn+1 



>N 



! ->N 



1 



>N 



where (j)' describe the new phases after the FD process. 
The diagram describes the situation just when the leader 
fires (first column) , the change in phases as a consequence 
of the emitted pulse (second) and finally the evolution 
of the system due to the linear driving up to the next 
firing (third). We believe that this is the simplest and 
most compact way to depict the process since we get rid 
of rotations that should be taken into account after the 
linear driving for any other relabeling method. Thus the 
transformation that describes this process reads 

$' =T n {$) = 1+M.J, 

where <fi' is a vector with N components since the zero-th 
component does not play any role in the description. In 
the above expression M n is an N x N matrix that can be 
written as 



(M 1 ) ii =5 i+1J -(l+e)S j , 1 



(2) 



J i+n,j 



+ eSj iSi+n,! Vn > 1. (3) 



In these expressions 6ij is the usual Kronecker delta. The 
sums should be interpreted modulus (N+l) and none of 
the subscripts can be either or N + 1. 

Since we are interested in emphasizing the mechanisms 
leading either to synchronization or pattern formation, 
we have considered very convenient to start our discus- 
sion with two illustrative situations where everything can 
be computed analytically and whose perfect understand- 
ing will help us to tackle the general case. 
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III. THREE OSCILLATORS 



0* =R A +M A -$*, 



(9) 



This is the simplest case which is worth analyzing, 
since the system formed by two units has been widely 
analyzed in the literature, see for instance |6| Jl^ , [l7|] . If 
we define a simple cycle as a sequence of firings in which 
each oscillator fires once and only once, there are only 
two possibilities for this system: 

• A: 0,1,2 

• B: 0,2,1 

provided oscillator is always at the threshold at the 
starting point of the dynamic evolution. Here, the num- 
bering corresponds to the firing sequence according to 
the initial spatial order in the lattice. Let us study both 
situations in detail. 



that is 



3 + 2e 



3 + 2e 



The stability of these fixed points is given by the ei 
values of which are 

3/2, 



-(H - 3)^ 2 ± i(n - l)fi^ 2 V4~]I 
2 

whose moduli are /i 3 / 2 . Depending on the sign of e 
fixed point is either stable (-) or unstable (+). 



(10) 

(11) 
gen- 

(12) 
the 



A. Order A: 0,1,2 

In this case, the sequence starts when oscillator fires, 
sending a pulse that changes the state of oscillator 1. 
Afterwards, the system is driven until oscillator 1 arrives 
to the threshold. According to our notation this process 
can be viewed as: 



1 



0^1 — n4>i = 02 

1 — > A^i— > 1 

2 — > 02 — ► 1 + 02 - M01 



(4) 



Then, we have transformed a state characterized by 
two phases 0i and 2 to a new one also characterized by 
two phases <f>' x and 2 such that (/)[ is always the phase 
of the oscillator that will receive the next pulse and keep 
this numeric order along the ring. In matrix notation the 
transformation can be written as follows 



-H 1 
-fi 



(5) 



The complete cycle is constructed by applying three 
times this transformation pi oT\oTi (0)] . In other words 



r = r a - 

The independent term Ra is 



Ra 



(6) 



(7) 



where 1 is a column vector of l's, 1 is the identity matrix, 
and the matrix of the transformation is defined as 



(8) 



From this expression it is easy to compute the fixed points 
of the transformation, which are solutions of the equation 



B. Order B: 0,2,1 

Now, oscillator 1 receives the pulse but it is oscillator 
2 the one that leads the driving and arrives first to the 
threshold 



1 -> -> 1 - 2 

4>\ — > ufa— > l + n<j>\ - 4>2 

02 — > 02 — > 1 



Therefore, the new phases are 



-1 
A 4 -1 



9i 

02 



(13) 



(14) 



Again the complete cycle is constructed by applying 
three times this transformation [Tg oT 2 o T2(0)]: 



0"' = R B +M b ■ 
where the independent term Rb now is 

R B = (l+M 2 + M 2 -M 2 ) -1 
and the matrix of the transformation 
M B = M 2 • M 2 • M 2 . 



(15) 



(16) 



(17) 



The fixed point of this transformation, which is solution 
of the equation 



— Rb 



us • 



is 



(18) 



(19) 
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tf = §±f (20) 

The stability of this fixed points is given by the eigenval- 
ues of Ms that are 

^-l±i(p-l)^T^T (2i) 

whose moduli are again /i 3 / 2 . Therefore, both fixed 
points describe the same physical behavior that is in- 
dependent on the particular order in which oscillators 
fire. 



C. Phase ordering 

In the previous subsections we have studied some fea- 
tures of the final state of the system when a sequence of 
identical transformations are applied successively. How- 
ever the reader can argue that those sequences are not 
the only possible dynamical evolution. Indeed, to deal 
with all the possible situations we should analyze what 
happens when a mixture of T\ and T2 are combined in an 
arbitrary manner to complete a cycle, and what sort of 
physical consequences derive from this fact. In addition, 
one may wonder whether an advancement can take place, 
i.e. if a given oscillator can fire twice before another ele- 
ment of the chain arrives to the threshold, breaking thus 
our definition of a simple cycle. Such issues are discussed 
in this subsection. To illustrate this point let us start by 
considering Fig. 1, where we have plot the evolution of 
the phases each time oscillator is at the threshold value. 
It is obvious that there are two different situations: a pos- 
itive or negative value of the coupling. In this figure we 
have analyzed the situation for e < 0. 

On the left hand side of Fig. I we have all the possible 
initial configurations and on its right hand side how they 
do transform when oscillator is again at the threshold. 
Regions A and B represent the sequences A and B de- 
scribed before, respectively. Here, we can see that states 
lying initially in one of these regions will approach the 
fixed points (attractors), since A and B are slowly shrink- 
ing. Therefore, once one starts with a given sequence no 
other one can be applied to describe the dynamical evo- 
lution of the system. The physical picture associated to 
the attractor fixed point is quite simple. The oscillators 
remain at a certain distance in the phase space (phase- 
locking). For larger dimensions (more oscillators) this 
fact induces the creation of complex spatial patterns. 




FIG. 1. Evolution of the phases <f>\ and cj>2, for a negative 
value of e. Right: initial configurations when oscillator is 
at the threshold value; left: new values of the phases when 
oscillator reaches the threshold again. The crosses corre- 
spond to the fixed points and each region corresponds to a 
given sequence of firings (see text). 

A special attention deserves region C which corre- 
sponds to a sequence of firings 010. This means that 
oscillator advances oscillator 2, a situation not covered 
before. Here we can see that the effect of this advance- 
ment is to reorganise the phases in such a way that after 
one cycle the old configurations fall in the basin of at- 
traction of region B and therefore, for them, sequence 
B must be applied forever. No more advancements can 
take place. The main conclusion is that advancements 
play a role only in the transient but not in the stationary 
properties of the dynamics. A more clear picture of the 
physical meaning of this fact will be provided in the next 
section. 

The mechanism for positive e is the opposite that for 
negative e. Regions A and B are enlarged every cycle, 
and the configurations move away the fixed points (re- 
pellers) until they cross some of the borders where at 
least two oscillators get absorbed and synchronize. Fig. 2 
shows this fact. The left hand side of the figure depicts 
the size of the basin of attraction of those configurations 
(closed lines) that when oscillator reaches the threshold 
again still require sequence A or B to evolve dynamically. 
The rest of the phase space is formed by states charac- 
terized by the fact that after the next firing of oscillator 
occurs at least two units will merge. These units are 
specified (underlined) on the right hand side of the fig- 
ure. Since two synchronized oscillators act as a single one 
and cannot be broken after a complete cycle, the problem 
now is equivalent to that of two oscillators. This dimen- 
sional reduction is the essence of synchronization. For 
mean-field models the word absorption has been coined 
to illustrate this phenomenon 0. 
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FIG. 2. The same as Fig. 1 for positive coupling. 

To address the question about the plausibility of hav- 
ing mixtures of consecutive transformations T\ and T 2 it 
is convenient to look at the problem from another per- 
spective. Instead of considering complete cycles it is bet- 
ter to analyze single firings. Let us suppose for simplic- 
ity negative (inhibitory) coupling between units. In this 
case, it is evident that the alternative application of both 
transformations lead to two possible options. The first 
case to be considered is the combination T\ o T 2 - We can 
observe that applying T 2 is inconsistent with the applica- 
tion of Ti afterwards; because the resulting configuration 
T 2 (4>) does not satisfy the possible different phase orders 
necessary to apply T\. On the other hand, the order 
T 2 o T\ implies one advancement between oscillators, cor- 
responding to region C (see Fig. 1) of phase space. This 
situation has been discussed previously to play the role of 
a transient dynamic behavior of the system. Thus we can 
conclude that, in general, the advancements will make 
the phases to be reordered until the system reaches a 
configuration which is consistent with only one sequence 
of transformations. 



IV. FOUR OSCILLATORS 

In this case there are 6 different orders for the oscilla- 
tors to complete a simple cycle: 

• A: 0,1,2,3 

• B: 0,1,3,2 

• C: 0,2,1,3 

• D: 0,2,3,1 

• E: 0,3,1,2 

• F: 0,3,2,1 

Analogously to the 3 oscillators case we define the ma- 
trices of the transformation between successive firings of 
two oscillators, according to the jump between oscillators 
that fire successively: 





M 3 




(22) 



(23) 



(24) 



We can easily compute the eigenvalues of those matri- 
ces. For Mi and M 3 the eigenvalues have moduli larger 
(smaller) than 1 for positive e (negative e). However, 
for M2 there is one eigenvalue with modulus equal to 1. 
This will be very important when discussing the stability 
of the fixed points. 

According to these jumps the transformation of a com- 
plete cycle for the different orders are constructed in the 
following ways 



• A: 0,1,2,3 

• B: 0,1,3,2 

• C: 0,2,1,3 

• D: 0,2,3,1 

• E: 0,3,1,2 

• F: 0,3,2,1 



Ti o Ti o Ti o Ti 
T 2 oT 3 oT 2 o Ti 
Ti o T 2 o T 3 o T 2 
T 3 o T 2 o Ti o T 2 
T 2 o Ti o T 2 o T 3 
T30T30T30 T 3 



For instance, let us write, in matrix form, the trans- 
formation for the case B: 



? = 1 + : 



1 + : 



ii 2 • - 



■ 1 + : 



12 • • . 



• 1+ 



+ M 2 • M 3 • M 2 • Mi • (f>. 



(25) 



We can now proceed to compute the fixed points asso- 
ciated to these transformations which are: 

• A: 1, 



4+3e' 4+3e' 4+3e ' 



• B- 1 i±£ 

±j. ±, 2+£ , u, 2+£ . 

• C: 2+1' l ' 2+7- 



• D: 1. 5Ti.l. 2+7- 



• E- 1 — i±l 



2+e' 2+£' 

F- 1 1 2+£ 3+6 
r ' ± ' 4+e' 4+£' 4+£- 
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Numerically, these fixed points are unique since this 
is ensured by the fact that all the eigenvalues of the 
matrix that multiplies <p arc different from 1 in all the 
cases. Nevertheless the fixed points for cases C and D 
are senseless since they do no verify the prescribed or- 
der, and hence they are physically unacceptable. Fur- 
thermore, the fixed points of cases B and C are physi- 
cally the same since they just differ in the order in which 
oscillators and 2 fire but they are synchronized. Thus, 
basically, we have to deal with fixed points where the 
phase difference between adjacent oscillators is roughly 
4>i+i — 4>i — 1/4, 2/4, 3/4(mod 1), with small corrections 
depending on e. 

Let us focus now on the process from the point of view 
of single firings again. Notice that Xi and T3 are trans- 
formations that can be applied alone successively leading 
to a natural complete cycle. However, T2 can not be ap- 
plied in this way because this would lead to an unphysical 
situation where only a group of oscillators will fire. This 
means that T 2 must be combined with Ti or T 3 in order 
to be physically acceptable, then the effect of having an 
eigenvalue with modulus 1, that can carry a metastabil- 
ity on the system, is avoided by this combination. The 
combinations of Ti with either T\ or T3, can be done in 
the way described by the orders B, C ,D or E, that always 
give rise to a chessboard type pattern. Other combina- 
tions can be obtained applying Ti o T 2 or T 3 o T 2 forever, 
this situation elicits advancements but nevertheless the 
resulting patterns are chessboard type again. 



V. N+l OSCILLATORS 

Although, in general, to construct a complete cycle is 
not a trivial mechanism we can infer some keys about the 
behavior of the system from the FD processes as we have 
done for the simplest lattices: compute the fixed points 
and their stability. 



A. Fixed points 

Due to their different behavior we will have to distin- 
guish again the cases n = 1 and n > 1. Thus, for the 
first situation, we have to solve 



01 
02 



M0i + 02 
M0i + 03 



i>jv-i = 1 - M0i 
i»jv = 1 - M0i 



9/V 



(26) 



Simply summing up all the equations we are left with 
4>i = N(l — ufa) and then we get 



N 



N + l + Ne' 



(27) 



Notice that the other phases at the fixed point can be 
obtained from this one, since 0j-i — j = 1 — /U0i V j = 
2, . . . , N. This fixed point corresponds to a situation in 
which all the oscillators fire in turn following their lattice 
ordering. 

On the other hand, for n > 1 we have to solve 



+ 0n+l 



9/Y 



t>N-n = 1 - 0r 
t>N-n+l = 1 — 0n 
t>N-n+2 = 1 - 4>n + M01 
f>N-n+3 = 1 - 0n + 02 



1 - , 



(28) 



Summing up again all the equations we now obtain 



i>i+0„=7V(l-0„)+M0i- 



(29) 



Now we notice that it is not enough to get one of the 
values of the phases. We have to close the system of 
equation by means of the following procedure: 

01 = 1 - 0„ + <?Wl = 2(1 - (j) n ) + 02n+l 

and so on. Again all the subscripts are understood mod- 
ulus N + l. This procedure is repeated until we reach 
4>N- n +i which closes the dependence of 4>\ on <p n . Ob- 
viously, a necessary condition to close it is that N + l 
and n do not have common factors. This procedure is 
iterated m n — 1 times, where m n verifies 

(1 + m n n) mod (N + 1) = 0, 

and it exists and is unique for each n < N+l. Then 
for a given N we will have to consider all the values of 
m n between 1 and N without common factors. We can 
therefore obtain that at the fixed point 



N ■ 



N+1 + m n e 



and 



N +1 + m n e 



(30) 



Again, this situation would correspond to a sequence 
of FD processes of jump n. It is easy to convince oneself 
that in both cases (n = 1 and n > 1) one can build a com- 
plete cycle by applying N+l times the transformation 
T n which does not change the fixed points. In principle, 
this successive application could make new fixed points 
to appear, but this is forbidden by the fact that the mod- 
uli of the eigenvalues of M„ is always larger than 1 for 
£ > 0, and smaller than 1 in the opposite case, whenever 
N+l and n do not have common factors, as we show in 
the appendix. 

But we still do not know what happens when N+l and 
n have common factors. As we show in the appendix in 
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this case there exists at least one eigenvalue of modulus 1 
and this fact does not ensure us about the existence of a 
solution for the set of algebraic equations; even when this 
solution can exist it usually gives rise to unphysical sit- 
uations, as, for instance, to values of the phases that are 
either below zero or above one. We showed explicitely for 
the case of 4 oscillators that under these circumstances 
the transformation of an FD process of this kind has to 
be combined with other transformations with noncom- 
mon factors and that this leaded to advancement between 
oscillators and to the formation of a chessboard like pat- 
tern. This is also what happens in the general case. The 
matrix associated to this combination will have eigenval- 
ues with moduli different from 1 and will guarantee the 
existence of the fixed point. Let us assume that a given 
spatial structure with period p — (N + l)/n exists and 
then there are some oscillators which are synchronized 

(4> = 4> P = 4>2 P = ■ ■ ■ = 4>{n-x)p\ 4>i = 4>p+\ = • • • ; • ••)■ 

Hence we need only to consider each spatial period, since 
the transformations of jump p correspond to oscillators 
that fire exactly at the same time. Within each period 
there will be different possibilities for the magnitude of 
the jumps n p and then the fixed points will be character- 
ized by 

rrir, nm„ 

* = ns_ = "p _ 31 

P + m np e N+\ + nm np e 

By combining Eqs. ©, ©, and ©, we realize 
that for a given N + 1 there always exists a value of to 
(0 < to < N + 1) such that 

^ = N+l + me ^ 

will be the phase of oscillator 1 at a fixed point. We have 
used this fact to identify the fixed points of the dynamics 
in simulations of lattices of a few oscillators, as we will 
see later on. 

In principle one can still think on the possibility 
of other fixed points corresponding to combinations of 
transformations not described above. For instance, suc- 
cessive applications of two transformations of different 
values of n; but this case will necessarily involve advance- 
ments between the oscillators, which, as we have already 
discussed, are only important in the transient but not in 
the approach to the final state. This, of course, can make 
the transients to become quite large, as we have observed 
in the computer simulations, but the only final states for 
an inhibitory coupling are those described earlier in the 
text. It is also important to note that depending on the 
strength of the coupling, and on the number of oscilla- 
tors, there will be some fixed points, for an inhibitory 
coupling, that will not exist; those that verify 

N + l 

(33) 



e < 1 



For instance, for the three oscillators case this happens 
for to = 2 when e < —0.5 [see (10)] and then region A 
disappears and the fixed point corresponding to B is the 
only possible final state. 



B. Stability of the fixed points 

After having shown the existence of the fixed points 
for single FD processes and extending this calculation 
to complete cycles, one needs to compute their stability. 
Since the calculation of the eigenvalues of the matrices is 
a lengthy but straightforward procedure we have left it 
for an appendix. There we show that the results easily 
obtained for three and four oscillators also apply to the 
general case, i.e. for e < the fixed points are attractors, 
whereas in the opposite case they are repellers. On the 
one hand, the attractiveness of the fixed points enables 
the formation of spatio-temporal patterns of phase-locked 
oscillators. This works not only for the structures with 
different phases but also for the periodic ones. On the 
other hand, when these fixed points become repellers it 
makes that neighboring oscillators synchronize and from 
that time on they will act as a single unit; this absorption 
(or dimensional reduction in our language) is iterated un- 
til the whole system acts as a single unit which completes 
the mechanism of the synchronization of the lattice mod- 
els with very-short range interactions we have analyzed 
through the paper. 



C. Computer simulations 

In order to check the validity of our results we have 
made computer simulations on lattices of a few oscilla- 
tors. In Tables I and II we represent the percentage of the 
structures the system formed by N + 1 oscillators reaches 
as a stationary state for two different values of e as a func- 
tions of to which stands roughly for the phase difference 
between neighboring oscillators times N+l, see (|3^). 
There are several results in these simulations that de- 
serve further comments. For instance we can notice that 
the oscillators tend to keep the maximum phase differ- 
ence. The chessboard like structure has the largest basin 
of attraction when the population has an even number 
of oscillators; whereas for an odd number of oscillators 
there are two peaks with the largest phase differences. 
However these results depend slightly on the strength of 
the coupling since the maximum percentage appears for 
the maximum phase difference and the larger the peak 
the larger the phase difference. Thus we can understand 
the different behavior for the two different values of e. 
Another difference concerns the reduction and, eventu- 
ally, the disappearance of the basins of attraction of the 
fixed points that correspond to large values of to. This 
fact also affects the time the system needs to reach the 
stationary state; for instance, for smaller values of e not 
only the jumps along one cycle are smaller but there are 
also more attractive fixed points. On the other hand, 
we have corroborated that for excitatory couplings the 
only possible final state is synchronization, no matter 
how long the transient is. 
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TABLE I. Percentage of the final states the system formed 
by TV + 1 oscillators reaches, for e=-0.1. The first column 
stands for N + 1 and the first row for m in (|3^), which ap- 
proximately corresponds to the phase difference times TV + 1 
between consecutive oscillators, ft is averaged for 1000 ini- 
tial random configurations picking each phase from a uniform 
ditribution between and 1. 
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4 5 
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10 


3 


61.4 


38.6 
















4 


21.2 


71.3 


7.5 














- 




5.5 


56.5 


37.2 


0.8 












6 


1.2 


29.1 


58.0 


11.7 0.0 












7 


0.0 


11.3 


51.7 


34.2 2.5 


0.0 










8 


0.0 


3.8 


32.1 


50.6 13.2 


0.4 


0.0 








9 


0.0 


1.1 


16.2 


47.9 30.8 


3.9 


0.1 


0.0 






10 


0.0 


0.2 


6.4 


34.6 43.8 


14.3 


0.7 


0.0 


0.0 




12 


0.0 


0.0 


0.7 


9.2 34.7 


40.0 


14.0 


1.4 


0.0 


0.0 


15 


0.0 


0.0 


0.0 


0.5 5.9 


23.0 


38.8 


25.2 


6.2 


0.4 




TABLE 11. The same as 


Table I for e~- 
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51.9 


48.1 
















4 


17.5 


67.3 


15.1 
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4.7 


46.7 


45.0 3.61 
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1.0 


22.8 


54.. 


3 20.7 


0.7 










7 


0.2 


8.4 


43.0 40.9 


7.4 


0.1 








8 


0.0 


2.8 


24.3 48.3 


22.5 


2.06 


0.0 






9 


0.0 


0.8 


11.35 39.7 


37.7 


9.7 


0.6 


0.0 




10 


0.0 


0.2 


4.4 26.0 


42.9 


23.1 


3.4 


0.1 


0.0 



VI. CONCLUSIONS 

In order to analyze the mechanisms of synchroniza- 
tion and the formation of spatio-temporal structures we 
have introduced a very simple model of pulse-coupled 
oscillators: a one-dimensional ring with unidirectional 
coupling. Despite this apparent simplicity it conserves 
all the features of low-dimensional systems subjected to 
short-range interactions which develop large-scale struc- 
tures. 

Although the dynamic evolution of the system involves 
two time scales, a slow one for the driving and a fast one 
for the interaction, we have constructed return maps that 
gives a complete information of the system. Concerning 
the maps, we have been able to compute exactly the fixed 
points of the dynamical evolution and their stability. 

For a negative (inhibitory) coupling the fixed points 
are attractors of the dynamics. Each one of these at- 
tractors has a well defined basin of attraction, although 
in some cases those regions are not simply connected. 
Since the evolution is discrete there are jumps among 
non-connected regions that correspond to advancements 
between oscillators. The advancements are only impor- 
tant in the transient dynamics, until the phases of the 
oscillators lie in the final basin of attraction of the fixed 
point. The volumes of the basins of attraction are differ- 
ent and depend on the value of the coupling, as we have 
checked by means of computer simulations of a few os- 
cillators. However, the states with the maximum phase 
difference between neighboring oscillators seems to be the 
preferred ones. 

On the other hand, for a positive (excitatory) coupling 
the fixed points arc rcpcllcrs of the dynamical evolution. 
Although in a configuration space with a multiplicity of 
repellers one can think that the system will jump from 
one region to another this is not our case. There are ab- 
sorbing barriers surrounding the repellers; when the sys- 
tem reaches one of these barriers it means that at least 
two neighboring oscillators have synchronized. When 
this happens the set of synchronized oscillators acts as 
a single unit that cannot be broken. From that time on 
we only need to consider a reduced number of units; we 
call this fact dimensional reduction since the new system 
can be described in terms of matrices with less compo- 
nents. This process of absorption is iterated until the 
system reaches a completely synchronized configuration. 

The present work only concerns the qualitative behav- 
ior of a population of pulse-coupled oscillators; never- 
theless, a quantitative behavior about the time a given 
population needs to reach the stationary state, either a 
synchronized one or a spatio-temporal pattern, would be 
desirable in order to complete the description of such 
systems. Another interesting question is related to the 
stability of the different structures with respect to fluc- 
tuations. 
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APPENDIX: 

In this appendix we will compute the bounds of the 
eigenvalues of the matrices M n defined in Eqs. (2)- (3). 
All the information will be extracted form the charac- 
teristic polynomial Pjy, n (X,e) which corresponds to the 

determinant |A1 — M„|. In this case we will distinguish 
several situations. First of all when N + 1 and n have 
common factors it is easy to see in the determinant that 
always exists a minor that corresponds to eigenvalues of 
modulus 1; anyway as we explain in the text we do not 
need to care about this case. 

The case n = 1 is very simple to compute. By simple 
inspection of the determinant it is easy to see that 

P Na (\,e)=\P N - 1A (\,e) + n (Al) 

and hence 

P N ,i(X, e) = X N + fiX"- 1 + fi\ N - 2 + . . . + fx\ + /i. 

(A2) 

In order to compute the case n > 1 it is convenient to 
introduce the following similarity transformation matrix 

where again the subscripts are understood modulus N+l, 
such that when it is applied to M n , M' n = ■ M n ■ B, 
one gets 

With this procedure we have converted a, in principle, 
complex matrix in a much simpler one. The new matrix 
has —1 in the last column and l's just below the main 
diagonal, unless at column m n and row m„ + 1 where it 
has 1 + e. In order to compute the eigenvalues we should 
notice that when e = the characteristic polynomial is 



P^ n (A,e = 0) = A Ar + A Ar - 1 +A Ar - 2 



+ A + 1. 



Then to compute it for e =/= we expand the determi- 
nant around the previous case and one realizes that the 
columns that are at the right and the rows that are below 
this element will not contribute. Thus one has 



N. 



t (A,£) 



X N + x N -' 



x N - 2 + . 



+ X m n + 



1) 



(A3) 



The last expression can be taken as general, i.e. including 
the n — 1 case, bearing in mind that mi = N. 



It is obvious that to compute the roots of ( |A3| ) is an 
unnecessary task, since the only needed information con- 
cerning the stability are the bounds of the eigenvalues. 
In order to compute these bounds we will look at the 
properties of the characteristic polynomial. It can be 
rewritten as 



N 



m n — l 



Pv,„(A, £ ) = ^A fc + £ Xk = 

k=0 k=0 

X N+1 _ 1 A m„ _ 1 



A — 1 



A - 1 



(A4) 



Then the eigenvalues of the matrix will correspond to the 
roots of 



(A 



N+l 



1) + e(A m " - 1) = 



unless the artificially introduced A = 1. These roots will 
verify 



N+l £X m 



1 + e. 



This modulus has upper and lower bounds given by the 
sum and the difference, i.e. 



|A| 



Af+l 



< l+s< |A| iV+1 + |e||A 



Thus for e > 0, we use 



lN+1 



e|A| m " > 1 + e 



which implies that |A| > 1. On the other hand, for e < 
we take 

|A| W+1 - |e||Ap" <l-|e| 

which, in turn, implies that |A| < 1. The final point 
is to show that the equality can only be fulfilled when 
both complex numbers have the same direction. It is 
easy to see that this can only happen when N + l and 
m n have common factors. Since this fact is avoided in 
this demostration we are left with the fact that for an 
excitatory coupling (e > 0) the eigenvalues are larger 
than 1 and the opposite for an inhibitory coupling (e < 
0), as we wanted to show. 
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